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We introduce a stochastic partial differential equation capable of reproducing the main features 
of spatiotemporal intermittency (STI). Additionally the model displays a noise induced transition 
from laminarity to the STI regime. We show by numerical simulations and a mean-field analysis 
that for high noise intensities the system globally evolves to a uniform absorbing phase, while 
for noise intensities below a critical value spatiotemporal intermittence dominates. A quantitative 
computation of the loci of this transition in the relevant parameter space is presented. 
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Spatiotemporal chaos (STC) is a complex behavior, 
common to many spatially extended nonlinear dynamical 
systems This behavior is characterized by a com- 

bination of chaotic time evolution and spatial incoher- 
ence made evident by correlations decaying both in space 
and time. In spite of considerable theoretical and exper- 
imental effort devoted to give a precise definition of STC 
and its different regimes, the present status is still un- 
satisfactory. A possible strategy to make progress in the 
understanding of STC is to investigate scenarios based 
on simple models, with few controlled ingredients, that 
reproduce the spatio-temporal structures under study. 
Such models are instrumental in searching for generic 
mechanisms leading to such complex behavior. Among 
the available scenarios, few of them consider the frame- 
work of stochastic partial differential equations (SPDE) 
to describe STC. A successful example is, however, 
the mapping of the Kuramoto-Shivashinsky equation (de- 
scribing a STC regime named phase turbulence) to the 
stochastic model of surface growth known as Kardar- 
Parisi-Zhang equation 

A particular instance of STC is a regime called spa- 
tiotemporal intermittency (STI) which is present in a 
large variety of systems Generally speaking, this 

regime is a chaotic spatiotemporal evolution (the turbu- 
lent phase) irregularly and continuously interrupted by 
the spontaneous formation of domains with a wide range 
of sizes and lifetimes, where the behavior is ordered (lam- 
inar) . The borders of the laminar regions propagate as 
fronts and eventually cause the collapse of the corre- 
sponding region into the turbulent background. There 
are strong indications that the STI regime has many fea- 
tures in common with phenomena of probabilistic nature. 
For example, it appears in some systems that critical ex- 
ponents at the onset of STI are in the universality class 
of directed percolation 13 1. STI has also been related to 
nucleation Q , another process associated with stochas- 
tic fluctuations. However, no description of STI in terms 
of SPDEs has been put forward so far. 

The purpose of this Letter is to introduce a model, en- 



tirely based on a simple SPDE, that describes the main 
features of STI, and report the existence of a noise in- 
duced transition from laminarity to STI. The laminar 
phase is associated to an equilibrium state called absorb- 
ing in the SPDE parlance. The role of the turbulent 
phase is played by a strongly fluctuating state driven by 
noise. For either large values of the noise intensity or 
small values of the diffusion constant, the system is glob- 
ally attracted to the absorbing state. However, for a fixed 
diffusion rate, STI sets in below a critical noise intensity. 
We investigate the nature of this transition both ana- 
lytically using a mean-field scheme, and numerically by 
means of indicators such as the change of the mean ve- 
locity of fronts between absorbing and turbulent regions, 
the one-site probability distribution function and the or- 
der parameter characterizing the average system fraction 
in the laminar phase. 

We consider the following one-dimensional Ito-SPDE 
fuj for the evolution of a real field u(x, t) fig]: 
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The r.h.s. of (|l|) is composed of a gradient term de- 
rived from a potential 4>(u), a diffusive spatial coupling 
with diffusion constant D, and a multiplicative Gaus- 
sian white noise £,(x,t) of zero mean and correlations 
(£(x,t)£(x',t')) =2 6(x-x')6(t-t'). For the results de- 
scribed below, the essential aspects of our model are (i) a 
bistable potential <fi(u) and (ii) a multiplicative noise func- 
tion g(u) which vanishes at the metastable state. Differ- 
ent forms of the noise amplitude g(u) give rise to different 
universality classes M . 

We make the simplest choice for the function g(u) = u 
Jl6| so that the sign of u(x, t) is preserved during the evo- 
lution of our system, allowing us to restrict our attention 
to the case u(x, t) > 0, by picking positive defined initial 
conditions. With this phase-space restriction, we choose 
the potential <fi(u) = au 2 — u A + hu 6 (with h > 0), which 
for < a < clm — 1 /4ft- has a relative minimum at u = 
and an absolute minimum at u = u+ separated by a max- 
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imum at u — u_ flT7|] , At a = <2m we are at the "Maxwell 
point" where <p(0) = tj>(u+), while at a = the u = 
state switches from metastable to unstable. In the ab- 
sence of noise a front between a region of u — and a 
region of it = u + moves at finite speed towards the for- 
mer, due to diffusive coupling. Consequently, any finite 
region of u = surrounded by u + will eventually shrink 
letting the system evolve into the uniform w + attractor. 
This is the case except for initial configurations laying 
completely on the relatively small basin of attraction of 
u = 0. 




FIG. 1. Space-time evolution of u(x, t) in a situation where 
the STI phase persists for all times (black: u — 0, grey: 
u > 0) . Space ranges in x £ (0, 400) (with periodic boundary 
conditions) and time in t £ (0,90). The numerical method 
uses a finite difference approximation for the Laplacian, and 
an Euler algorithm for the time integrator which respects the 
positivity of the field. The initial condition is chosen ran- 
domly in the interval uo(x) £ (0,2.4). Other values of the 
parameters are e = 0.95, a = 0.5, D = 2.0, h — 0.22, Ax = 1, 
At = 0.001. 

In the presence of multiplicative noise, the state u(x) 
becomes an absorbing state |?]], i.e. a state where the 
system can be driven to by fluctuations but not removed 
from by neither these nor by the deterministic dynam- 
ics. This is because the selected noise amplitude func- 
tion g(u) vanishes at the fixed point u — 0. On the other 
hand, when the field takes values near u_|_, the absolute 
minimum of 4>{u), fluctuations are always active. How- 
ever, without diffusion (D = 0) the fate of the system is 
to asymptotically settle in the absorbing state, since for 
any e/Oa fluctuation large enough to push the system 
to the basin of attraction of u = will eventually occur. 

From this analysis, we infer that noise and diffusion 
play opposite roles on the asymptotic space-time evolu- 
tion of the field u. While noise nurtures the development 
of regions dominated by the absorbing state, diffusion 
favors the dominance of u values fluctuating around the 
global minimum of 4>(u). We associate the quiescent uni- 
form absorbing state of this system with laminarity and 
the disordered fluctuations around the global potential 
minimum with turbulence. In Fig. [l] we show a space- 
time plot of a numerical solution of Eq. ([!]). Notice that 



for these values of e and D there are occasional nucle- 
ations of regions of the laminar or absorbing state that 
eventually collapse under the progress of the fluctuating 
or turbulent state. The wide range of sizes and lifetimes 
of the laminar areas, characteristic of STI, is evident in 
the picture. 

A similar mechanism, also involving the nucleation of a 
laminar state, has been invoked to be responsible of STI 
in a deterministic reaction-diffusion system (Oj. In this 
case, the role of noise is played by a chaotic dynamics 
generated via a Hopf bifurcation. 
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FIG. 2. Space-time plot of u(x,t) starting from a random 
initial condition uo(x) of mean u+ for x £ (0,400) except 
in the interval x £ (150,250) where uo(x) — 0. The noise 
intensity e is changed during the evolution: e = 0.005 for 
t £ (0,30), e = 0.145 for t £ (30,105) and e = 0.250 for 
t £ (105, 180). Other parameters are a = 1.0, D = 2.0. 

In Fig. ^|, for fixed a and D, we compare the evolu- 
tion of fronts separating absorbing and active regions for 
three different noise levels. For the lowest noise intensity 
e these fronts invade the absorbing phase, as in the deter- 
ministic system, and the system asymptotically tends to 
be globally in the active phase u+. On the other extreme, 
for the highest e, these fronts invade on average the ac- 
tive phase which consequently tends to disappear. The 
evolution for an intermediate value of e shows a situation 
where the fronts are at rest on average This sug- 

gests a possible mechanism for a noise induced transition 
from the STI phase to the absorbing phase. 

In Fig. ||a we give a quantitative evidence of this tran- 
sition. There the average fraction R of the system in 
the absorbing phase is plotted as a function of the noise 
intensity e. This figure shows a transition from R w 
(STI) to R = 1 (absorbing phase) at roughly the same 
critical noise intensity as for the front velocity reversal. 
Another useful measure to characterize this transition is 
the one-site stationary probability density function (pdf) 
0, P{u). In Fig. IJd we plot P(u) for different values of 
the noise intensity. We observe a transition from a pdf 
which has a hump around the active phase u+, to a pdf 
which is highly peaked at the absorbing state u — 0. No- 
tice that this transition occurs at the value of e for which 
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lal ( a = - 5 ) 



P(u;u) 



1.0 - 



°< 0.5 



0.0 







J 





0.80 0.90 1.00 1.10 




0.10 




FIG. 3. (a) Fraction portion of absorbing phase R as a 
function of e and two different values of a. The data are the 
result of averaging over 40 realizations, for a time t = 900. 
Other simulation parameters as in Fig. hi (b) One-site proba- 
bility distribution P(u) for different values of e. The ordinate 
axis of the pdf corresponding to e = 1.025 is scaled down 10 
times. 

The transition from the active phase to the absorb- 
ing phase that we have numerically characterized above, 
can be described within a Weiss-like mean field the- 
ory specially devised to study the effect of fluctuations 
|jl9|| . Wc hrst notice that, by a suitable rescaling of the 
space variable x, (Q) can be rewriting as: Ut — f(u) + 
D u xx +g(u) £ with f{u) = -§£ and D = D/e 2 . We 
next consider a spatial discretization of the field where 
Ui = u(xi,t). One can write the multi-variate Fokker- 
Planck equation for the set of variables {u\, U2 ■ ■ .} which, 
after integration of all the u, variables except one, yields 
the equation for the stationary one-site pdf P(u): 

d u [(-/(«) + 2D (u - E(u)))P] + dl(g 2 (u)P) = (2) 

where E(u) = J u' P(u'\u) du' is the steady-state con- 
ditional average of the field at a nearest neighbor site of 
a site in which the field takes the value u. The above 
mentioned mean-field approximation takes E(u) ~ u, 
the yet unknown mean field value of u. This is analo- 
gous to the traditional Weiss mean-field approach in the 
theory of critical phenomena. We integrate (0) to obtain 
the steady state pdf: 



Z g{u) 



exp 



u f(v) - 2D (v - u) 
g(v) 2 



dv (3) 



with Z the normalization constant and where the depen- 
dence on the mean value u has been explicited. The value 
of u arises from the consistency relation: 



v P(v; u) dv 



(4) 



We can easily solve (Q) in the no-coupling limit Dq = 0, 
where P(u) ~ u ~ 2 ( 1 + a ) which is non-normalizable (re- 
member that a > 0) so in fact we have P(u) = S(u), and 
hence u = 0, which describes the absorbing phase. On 
the other hand, the limit Do = oo can be treated us- 
ing a saddle-point expansion in Dq which results in the 
equation f(u) — P0|| . This equation coincides with the 
steady state result of the deterministic analysis for a spa- 
tially averaged field. Such analysis predicts a transition 
from the absorbing state u = to the state u = u + at 
the Maxwell point a = clm- These limiting results show 
that the transition from STI to absorbing appears in our 
parameter regime < a < bai as a joint effect of fluctua- 
tions and spatial coupling. We note that other transition 
to an absorbing phase, recently studied in the context of 
SPDE's, already exists for D Q = @. 
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FIG. 4. Phase diagram of the noise induced transition from 
STI to absorbing phase. The symbols correspond to numerical 
simulations of (hi) for different values of D and e. The solid 
line comes from the solution of the mean-field consistency (^) 
[§• 

In order to obtain results for intermediate values in the 
coupling Do we solve numerically the consistency relation 
([|). In Fig. U we show a phase diagram where the mean 
field approach and the numerically determined transition 
from (^), are compared. 

A possible variation of our model is to include a finite 
correlation length in the Gaussian white noise (Fig. ||). 
This modified model recreates the characteristic triangu- 
lar structures displayed by several two-component deter- 
ministic reaction-diffusion equations [ pd|p2| . Note that 
in this case the transition to the laminar state occurs 
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for non-zero velocity of the fronts, in analogy to what is 
observed in deterministic models ItOJ. 




(£(x, t)£(x' , t')) = exp(— — x I )S(t — if) for two noise intensi- 
ties: e = 0.35 for t G (0, 150 7 ) and e = 0.50 for t G (150, 250); 
x G (0, 400) and 7 = 8.0, D = 2.0. 



As a conclusion, we have shown that STI can be mod- 
eled by means of a simple SPDE. Our model provides a 
new insight into the complex behavior of the STI phase as 
a stochastic nucleation of an absorbing metastable state. 
One of our main results is the occurrence of a noise in- 
duced transition where STI disappears in favor of an ab- 
sorbing phase at sufficiently high noise intensity or low 
enough spatial coupling. Notice that both, transitions 
to absorbing states Hj7j and noise induced phase tran- 
sition |2^] are known in the SPDE context. The former 
phenomenon does not require spatial coupling. The lat- 
ter has been so far associated with transitions between 
statistically homogeneous stationary states. The noise 
induced transition shown here is a genuine consequence 
of stochasticity in extended systems. It describes a tran- 
sition between an absorbing phase and a dynamically ac- 
tive and structured phase which is a form of STC. Specific 
features of this active phase are described by considering 
space-time configurations associated with the individual 
realizations of the stochastic dynamics. In summary we 
believe the SPDE approach to be general to explore new 
generic mechanisms of spatiotemporal chaos. 
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